1 



AI =1/2 rule from staggered fermions. 

D. Pekurovsky and G. Kilcup a 

department of Physics, the Ohio State University, 174 W. 18th Ave., Columbus OH 43210, USA 



00 

On 
Or 
in 



> 
in 



On 
O 
0O 

D 



X 



We present our latest results for the AI — 1/2 rule, obtained on quenched ensembles with /3 = 6.0 and 6.2, 
and a set of JV/ = 2 configurations with /3 = 5.7. The statistical noise is quite under control. We observe an 
enhancement of the A J =1/2 amplitude consistent with experiment, although the systematic errors are still large. 
We also present a non-perturbative determination of Zp, Zs and the strange quark mass. We briefly discuss our 
progress in calculating e'/e . 



1. Introduction and methods 

It is well-known that the AI =1/2 channel 
of non-leptonic kaon decays is enhanced com- 
pared to the AI — 3/2 one. In particular, 
lu = ReA /ReA 2 = 22, where A 0<2 e iS °- 2 = 
((tttt)i=o,2\Hw\K°) This talk is a status report 
on our work in calculating these matrix elements 
(MEs) using staggered fermions. We have also 
computed MEs relevant for e'/e . 

The effective weak Hamiltonian for this prob- 
lem is as follows: 



H eff 



-j=V ud V* s 



10 

^2\Zi(jt) + Tyi(lJ,) 0;(/i),(l 
i=l 



where r = -V td V t */V ud V* s , Z{ and are Wilson 
coefficients. We work in the standard basis of the 
10 four-fermion operators defined in 

Our task is to compute (mr\Oi\K). Putting 
a two-pion state on the lattice is a well-known 
technical problem. We use the chiral perturba- 
tion theory method || to relate (mr\Hw \K) to 
(tt\H w \K) and (0\H w \K). This is equivalent to 
subtraction of a single lower-dimension operator 
O su b = {m s + m d )sd + (m d ~ r7i s )s7 5 d, which is 
the only operator allowed to mix due to the chiral 
properties of staggered fermions. This procedure 
does not take into account the higher order cor- 
rections in chiral perturbation theory (ChPT), in 
particular the final state interactions of the pi- 
ons. These corrections are known to be large, 
which introduces a big systematic uncertainty in 
our results. 



We follow the strategy of ME computation 
with staggered fermions yj] and compute three 
types of fermion contractions, known as "eight", 
"eye" and "annihilation" diagrams. The latter 
two types are quite noisy, but we gained enough 
statistics to bring the noise down to an acceptable 
level for all basic operators 0\-0\q. 

Table 1 shows our simulation parameters. Our 
lattice is replicated 4 times in time direction to 
avoid contamination from excited states. We use 
degenerate mesons (m s = m d — m„) and gauge- 
invariant, tadpole- improved operators with stag- 
gered fermions. 
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Figure 1. Bk vs. meson mass squared. Vertical 
line here and in other plots indicates the physical 
kaon mass. The fit is of the form y = a + bx + 
cx logx. 
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Table 1 

Simulation parameters 



N f 





Size 


L, fm 


N 





6.0 


16 3 x (32 x 4) 


1.6 


216 





6.0 


32 3 x (64 x 2) 


3.2 


26 





6.2 


24 3 x (48 x 4) 


1.7 


26 


2 


5.7 


16 3 x (32 x 4) 


1.6 


83 



2. ReA 2 and (O k ) 

ReA 2 can be related by ChPT to (K°\O K \K ), 
so at the lowest order we just need to com- 
pute f K and B K (denned in (K°\0 K \K°) = 
8/3 m 2 K f K B K , where K = sj L d sj L d). 
The Bk parameter involves calculating only 
"eight" contractions and is now well studied 
(e.g. [||J|]). The form of the chiral be- 
haviour Bk = a + bm 2 K + c m K log m K || pro- 
duces a reasonable fit and gives a finite non- 
zero value in the chiral limit (Fig. 1). However, 
the physical ReA 2 (proportional to (Ok)/iti k = 
8/3 J k Bk) is very sensitive to the meson mass, 
contrary to a naive expectation (see Fig. 2). 
This is due to a sizeable slope in Jk vs. m 2 K 
(which is of the same order as the physical 
slope) (Fig. 3). Thus there is a large un- 
certainty due to unknown higher order terms 
in ChPT. Naively taking the meson mass 
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Figure 2. KeA 2 vs. meson mass squared. The fit 
is of the form y = a + bx + cx log x. The vertical 
line indicates the physical kaon mass. 
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Figure 3. /k vs. meson mass squared. 

M 2 = (m 2 K + ml )/2 and using our quenched 
value of Bk in continuum limit yields: ReA2 = 
(1.67 ± 0.02) • 10" 8 GeV, to be compared with 
experimental 1.25 • 10 -8 GeV. 

3. ReAo and AI = 1/2 rule 

Results for Re^4o f° r quenched (3 = 6.0 and 
6.2 ensembles are shown in Fig. 4. The depen- 
dence on the meson mass is much smaller than 
that of KeA 2 . We have checked the lattice vo- 
lume dependence and found it to be small for lat- 
tice sizes 1.6 fm and above. However, the results 
significantly depend on the lattice spacing (see 
Fig. 5). The (3 = 6.2 point seems to bring the 
continuum value below the experiment. However, 
final state interactions could additionally raise it 
by as much as 100%. Finally, we checked the 
effect of unquenching and found it to be small 
compared to noise (see Fig. 6). 

We show the ratio Re^b/ReAo for quenched 
{(3 — 6.0) and dynamical data sets in Fig. 7. The 
result is of the same order of magnitude as expe- 
riment. However, the dependence on the meson 
mass is so large that it prevents us from a more 
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Figure 4. Re^4o y S- meson mass squared for 
quenched (3 — 6.0 and 6.2. Finite volume study 
for j3 = 6.0 is also shown. 
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Figure 5. ReAo vs. lattice spacing squared. The 
cross is the result from the larger volume f3 = 6.0 
ensemble. 

definite conclusion. Clearly, we need higher order 
ChPT terms to elucidate this subject further. 

4. Operator matching and e'/e perspective 

Lattice and continuum operators have to be 
matched, for example by perturbation theory: 
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Figure 6. The effect of unquenching on Re^o- 
The lattice spacings of quenched and dynamical 
ensembles are comparable. 



One-loop perturbative matching works well for 
operators 0\ and O2 (relevant for Re^o): the 
corrections are small, and so is the q* depen- 
dence. However, for operators O5 - Os (relevant 
for e'/e ) the situation is much worse: several 
of the perturbative coefficients have not yet been 
calculated, while others (notably Cpp) are too 
large at one-loop order to trust the perturbation 
theory Thus, a non-perturbative matching pro- 
cedure is necessary to calculate e'/e . 

As a first step towards this procedure we have 
computed bilinear renormalizaton factors Z$ and 
Zp using a non-perturbative method based on the 
strategy by Martinelli et al. (see Fig. 8). Our 
best value for the strange quark mass in conti- 
nuum limit obtained with non-perturbative Zg 
in MS at 2 GeV is 103 ± 8 MeV, which is quite 
close to the one- loop result. 

Zp can be used to get a reasonable estimate of 
the renormalization of operators Oq and Os- The 
value of (Oq) obtained in this way is very different 
from the tree-level value: it is much smaller, close 
to zero. This would produce a negative e'/e , 
contrary to experiment. However, to give a more 
definite prediction we need to perform a full non- 
perturbative renormalization procedure. 



Or*(«*)= 0« + ^gjl^( 7i ,ln(^) 5. Summary 

•? We have obtained a reasonable statistical pre- 

+C ij )Of t + 0(g 4 ) + 0(a n ) (2) cision in studying ReA and ReA 2 as well as all 
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Figure 7. Isospin amplitude ratio vs. meson mass 
squared for quenched and dynamical ensembles 
with comparable lattice spacings. This enormous 
dependence on the meson mass comes entirely 
from the behavior of Re^2 (see Fig. 2). 
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Figure 8. Non-perturbative renormalization fac- 
tors Zs and Zp in MS vs. momentum scale /i for 
f3 = 6.2. The curve is the one-loop perturbation 
theory expectation (which does not distinguish 
between Zs and Zp in the chiral limit). 



matrix elements needed for s' /s . The biggest un- 
certainties in our prediction of the Al =1/2 rule 
are higher order ChPT terms and continuum ex- 
trapolation. In addition, a non-perturbative ope- 
rator matching needs to be done for e'/e . 
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